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Recently the dynamically assisted Schwinger mechanism, i.e., electron-positron pair production from vacuum by a combination 
of laser pulses with different time scales has been proposed. The corresponding results, which suggest that the rate of produced 
pairs is significantly enhanced by dynamical effects, are verified. Employing the framework of quantum kinetic theory intrinsically 
enables us to additionally provide momentum space information on the generated positron spectrum. 
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1. Introduction 

The description of a charged particle within Quantum Elec- 
trodynamics (QED) necessitates to take the corresponding an- 
tiparticle into account. Historically, the first example has been 
the prediction of the positron, the antiparticle of the electron. 
One may exploit the similarity of quantum field theory to many- 
body physics by introducing the Dirac sea picture: The vacuum 
is modeled by a sea of fermions with all negative-energy levels 
filled. The antiparticle is then represented by a hole in this sea. 
This makes plain that an antiparticle-particle pair can be created 
as soon as enough energy is available to bridge the mass gap. 

It became evident soon that one can create an antiparticle- 
particle pair either perturbatively from photons in high-energy 
reactions or non-perturbatively in the presence of very strong 
electric fields. The latter type has been called Schwinger mech- 
anism [1] despite the fact that it has been discussed already 
many years before 2,3]. This effect has attracted a lot of inter- 
est recently fl B UWM H OH [H d Q since laser facilities 



X 



such as the Extreme Light Infrastructure (ELI) or the European 
XFEL might possibly reach the required field strengths. 

Being aware of the possibilities of perturbative versus non- 
perturbative production mechanisms, the question arises how to 
switch in a quasi-continuous way from one to the other. In this 
respect it turns out that time-dependent strong electric fields do 
this when varying the relevant frequency. Consequently, allow- 
ing for different and suitably chosen frequencies leads to a kind 
of dynamically assisted Schwinger mechanism 

The aim of the investigation reported in this letter is twofold: 
First, the mechanism as such is verified employing a quantum 
kinetic equation instead of the (semi-classical) world-line for- 
malism. Second, and in view of the experiments even more im- 
portant, momentum spectra of the produced pairs are predicted. 
Whilst the considered case is certainly not realistic in the sense 
of representing experimentally relevant superpositions of elec- 
tric fields, it already demonstrates very nicely what kind of vis- 
ible effects in the momentum spectra are to be expected in ap- 
propriately tuned upcoming experiments. 



This letter is structured as follows: In Sect.|2]we set the theo- 
retical basis for our investigations and introduce a simple model 
for dynamically assisted Schwinger mechanismQ which is based 
on the proposal [7]. In Sect. [3]we discuss our numerical results 
on the momentum spectrum as well as on the particle density. 
We show that combining the mechanisms of perturbative and 
non-perturbative pair creation results in strong non-linear be- 
havior which in turn leads to an enhancement in the particle 
yield. 

2. Pair Production in Multiple Time Scale Fields 

2.1. Homogeneous Background Field Approximation 

With the advent of high-intensity laser facilities such as ELI 
or the European XFEL, we hope that it will become possible to 
create field strengths of the order of S c = m 2 /e ~ 10 15 V/cm in 
the focus of colliding laser pulses. A laser should be described 
as an ensemble of photons in principle, showing both temporal 
and spatial variation. However, due to the fact that a complete 
description of this system is not yet feasible, we will make sev- 
eral reasonable simplifications in our calculations: 

• We consider very high field strengths. Consequently, it is a 
reasonable assumption to treat the laser as classical back- 
ground field. We thus neglect processes which could lead 
to a possible depletion of the laser due to electromagnetic 
cascades as indicated recently 

• The Schwinger mechanism is driven by the electric field. 
Having in mind a situation of colliding laser pulses with 
canceling magnetic fields, we will ignore the magnetic 
field completely; for studies including magnetic field 
effects see e.g. 1 18, 1^1 . 



For a detailed and self-contained description, the interested reader may 
refer to 0. 
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The electric field is space- and time-dependent in general. 
Nevertheless, since the achievable spatial focusing scale is 
orders of magnitude larger than the Compton wavelength, 
we will ignore any spatial dependence. For studies on the 
spatial dependence of the Schwinger mechanism see e.g. 



Consequently, we assume a spatially homogeneous electric 
field which is represented by a vector potential in temporal 
gauge A^(t) = (0, A(f))- The spatial part is chosen to point into 
the e3-direction: A(f) = (0, 0, A(f)), such that the magnetic field 
vanishes, B(f) = V x A(f) = 0, and the electric field reads: 



E(r) 



dA(t) 
dt 



(0,0, £(f))- 



(1) 



2.2. Quantum Kinetic Formalism 



We employ the quantum kinetic formalism in order to obtain 
not only production rates but also to gain information on the 
momentum spectrum of created particles. This additional in- 
formation is highly relevant, particularly with regard to a possi- 
ble future detection of the Schwinger mechanism. The spectral 
information is encoded in the distribution function /(q, t): It 
has to be stressed that there is no clear interpretation of /(q, f) 
as long as the electric field is present. At finite times it can- 
not be considered as distribution function of real particles but 
only as a mixture of real and virtual excitations. Consequently, 
/(q, t) might be interpreted as the momentum distribution for 
real particles only at asymptotic times t — > +oo, when the ex- 
ternal field is switched off. This interpretation is also supported 
by S-matrix theory which clearly states that the identification of 
excitations of quantum fields with its corresponding particles is 
only possible at asymptotic but not at intermediate times. Alter- 
native interpretations, which do not account for this peculiarity 
of /(q, t), lead to somewhat astonishing results Bo, 12411 . 

The equation of motion for /(q, f) can be derived from canon- 
ical quantization, quantizing the spinor field fully but consider- 
ing the electromagnetic field as given background ISI2SI27I]. 
Alternatively, a derivation based on a Hartree approximation of 
the relativistic Wigner function for QED is possible I123[ 12811 . 
Due to the fact that we are mostly interested in the sub-critical 
field strength regime & <K £ c where the expected densities are 
rather low, we neglect collisions of created particles lfl9l 291. 



Additionally, the self-consistent field current due to created 
particles, which leads to an internal electric field and the de- 
crease of the background, should be taken into account in prin- 
ciple. Detailed analysis, however, showed that these contribu- 
tion can be safely neglected in the sub-critical field strength 
regime £ <§: & c as well |[l9, 30[ 31 ]. With these simplifications, 
the quantum kinetic equation for /(q, t) reads: 



/(q,?) = W(q,f) 



1 

j dt'W(q,t')[l 



•/(q,f')] cos [20(q;*,O] 



(2) 

where f(q, t) accounts for both spin directions due to the ab- 
sence of magnetic fields. Here, 



W(q, t) = eE( ' )£± an d 0( q; h , t 2 ) = f dt'<o(q, t') , 
w 2 (q,f) J 

h 

with q = (q ± , ^3) being the canonical momentum; e denotes the 
electron charge; e\ — m 2 + q\ is the transverse energy squared, 
whereas to 2 (q, t) - e]_+ pi{t) 2 is the total energy squarecQ Al- 
ternatively, this equation may be expressed as linear, first order, 
ordinary differential equation (ODE) system [f30h : 

f(q,t) = W(q,t)v(q,t) 

v(q,0 = W(q,t)[l~f(q,t)]-2<o(q,t)u(q,t) (3) 
u(q, f) - 2w(q, f)v{q, f) 

with initial conditions /(q, f,) = v(q, f,) = u(q, f,) = at f,- — > 
-00. Both the integro-differential equation (f2]i and the ODE 
system (0) are valid for any time-dependent electric field. 

In addition to the spectral information, the second important 
observable is given by the asymptotic particle number density 
n[e + <?~], which is calculated directly from /(q, t): 



n[e + e ] 



d 3 q 



f(q, °°) 



(4) 



2.3. Keldysh Parameter 

In order to create electron-positron pairs out of the vacuum, 
it is necessary to overcome the threshold energy of 2mc 2 from 
the filled Dirac sea to the unoccupied energy levels. Basically, 
there are two mechanism to achieve this: 



Schwinger pair creation |l|,|3j]: Virtual electron-positron 
fluctuations may gain the necessary energy e&L > 2mc 2 
by acceleration over a distance L in a static or slowly vary- 
ing electric field &. This effect might be considered as 
tunneling effect and is as such exponentially suppressed. 



• Dynamical pair creation B32l 13311 : In a time-dependent 
electric field with characteristic frequency u, the threshold 
energy might be provided by a single photon Hcj > 2mc 2 . 
For ha> < 2mc 2 , this effect might still occur by absorbing 
a multiple number of photons, but then becomes strongly 
suppressed as well. 

An important parameter facilitating this characterization is the 
Keldysh parameter, first introduced in the context of atomic ion- 
ization 13411 but later also taken over in the context of electron- 
positron pair creation from vacuum [32]: 



r 



T 



m 

e&T 



(5) 



Here, tj = m/e&is the characteristic time of a tunneling event. 
On the other hand, t denotes the characteristic time scale of the 
applied electric field. The ability to tunnel is facilitated by the 
electric field & which tilts the potential around the electric field 
axis, as shown in Fig.Q] 



2 The well known connection between canonical and kinetic momentum is 
p(f) = q - eA(t). 
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Figure 1: Left: Relativistic dispersion relation for vanishing electric field. 
Right: The electric field leads to a tilt in the relativistic dispersion relation such 
that tunneling out of the filled Dirac sea to the continuum is possible. 




Figure 2: Superposition of two Sauter-type electric fields Eq. {8} for to = 0. 
The relevant scales in this system are the larger field strength £i and the shorter 
time scale T2 . 



Consequently, the Keldysh parameter discriminates between 
two regimes: 

• Adiabatic regime y <k 1 (jj <K r): Tunneling through 
the barrier is possible since the electric field might be 
considered as static at any instant of time. In this 
regime, electron-positron pairs are created spontaneously 
by means of a non-perturbative tunneling process. 

• Anti- adiabatic regime y » 1 (jj s> r): Tunneling through 
the barrier is not possible because the shape of the electric 
field changes much faster than the time needed to tunnel. 
In this perturbative regime, electron-positron pairs might 
only be created by absorption of the necessary threshold 
energy and not via a tunneling event. 

2.4. Electric Field with Two Time Scales 

Multi-photon pair creation has already been experimentally 
verified at the SLAC E-144 experiment H35I1 . whereas the 
Schwinger mechanism has not been observed so far due to its 
strong exponential suppression of the order of exp(-7r£ c /£). 
Hence, various mechanisms have been proposed in order to 
partly overcome this suppression of the Schwinger mechanism 
by taking advantage of pair creation in the multiphoton regime 

The idea of the dynamically assisted Schwinger mechanism 
is rather simple [7]: Superimposing a strong but slowly varying 
electric field with a weak but rapidly varying electric field leads 
to a decrease of the tunneling barrier and, consequently, to a 
drastic enhancement of the electron-positron pair creation rate. 
As a very simple example we consider two Sauter-type electric 
fields 



Ei(t) = Si sech z (f/T,-) . 



(7) 



3 Note, that due to the fact that the Dirac equation is exactly solvable in terms 
of hypergeometric functions in a single background field Eq. (7), we are able to 
find an analytic expression for the distribution function f(q, t) at both finite and 
asymptotic times [ 23]. For simplicity, we only give the result of the asymptotic 
momentum distribution here: 



f(q, oo) : 



2 sinh \n \e£r 2 + fi - vl) sinh (n \e&T 2 - /j + yj 
sinh(27T/j) sinh(27rv) 



(6) 



with 



^ yjm 2 + q 2 ± + (q 3 - e&r) , 
T - ^m 2 + q 2 ± + (q 3 + e£r) . 



In our calculations we take E\(f) as the strong field, whereas 
£2(0 is considered as the weak field. Consequently, we have 
£2 «c £1 and T2 <k Ti . In terms of the Keldysh parameter we 
have y\ <k 1 and 72 s> 1 . We then consider the superposition 
of two such electric fields (cf. Fig. [2J: 



E(t) = £1 sech 2 (f/i-i) + £2 sech 2 ((f - r )/r 2 ) 



(8) 



The parameter to accounts for a possible shift of the individual 
pulse maxima with respect to each other. 

It must be stressed that the relevant Keldysh parameter for 
this problem is neither y\ nor 72 but a combination of the dom- 
inant scales of the two different electric pulses [7]. These rel- 
evant scales are given by the larger field strength £1 and the 
shorter time scale T2: 



m 



y 



e&\T 2 



(9) 



This combination of different scales results in a subtle change in 
the pair production process by means of a non-trivial combina- 
tion of Schwinger and dynamical pair creation. Consequently, 
the resulting distribution function /(q, f) will not be given by 
a naive addition of f\ (q, f) for the Schwinger mechanism and 
/2(q, t) for dynamical pair creation, as we will see in the next 
section. In our calculations we consider a given strong electric 
field E\ (t) in the adiabatic regime and a weak electric field £2(0 
which can be varied by changing either the field strength £2 or 
the pulse length T2. 



f[q] 
0.00003 - 




Figure 3: Comparison of the distribution function /(q, 00) (solid line) with 
/lCq, 00 ) (dotted line). The positions of the maxima of the distribution func- 
tions are marked by m and m\, respectively. The parameters are given in the 
text. /2(q, 00) is not shown because it is orders of magnitude smaller. 
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3. Numerical Results 

We now turn to the numerical results for the dynamically as- 
sisted Schwinger mechanism in an electric field with multiple 
scales. In a first instance we focus on the asymptotic distribu- 
tion function /(q, oo) of created particles. Subsequently, we dis- 
cuss the enhancement in the number density n[e + e~] of created 
particles in a specific parameter region due to the non-trivial 
combination of Schwinger and dynamical pair creation. 

3.1. Momentum Spectrum 

In Fig. [3] we consider the dynamically assisted Schwinger 
mechanism for a given setup, with the parameters for the strong 
pulse {fii = 0.25£ c ,ti = 10~ 4 eV~'} and the weak pulse {£2 = 
0.025£ c ,r 2 = 7 • \0- 6 eV- 1 }. For these parameters, y x * 0.078 
is deep in the non-perturbative regime whereas j2 ~ 11.18 is 
well in the multi-photon region. In this figure, we compare the 
momentum spectrum for a single Schwinger-like pulse with the 
momentum spectrum for a combined pulse Eq. <j8j for fo = 0. 
We should first make two remarks: 

• The parameter sets {£i,ti} and {£2^2} are chosen such 
that the contribution due to Schwinger pair production 
fx (q, °°) is orders of magnitude larger than that due to dy- 
namical pair production ^(q, 00). In order to observe the 
dynamically assisted Schwinger mechanism, it is of great 
importance to consider a parameter region which is in this 
sense dominated non-perturbatively (Schwinger-like). 

• The distribution function for single pulses is easily found 
by the analytic formula Eq. (|6}. This formula also served 
as a numerics benchmark throughout this work. 

We observe, that the momentum spectrum for the com- 
bined pulse /(q, 00) is not given by naively adding f\ (q, 00) 
and /2(q, °o) but shows a strong non-linear behavior. In fact, 
/(q, 00) and f\ (q, 00) nearly coincide for small kinetic momenta 
but differ for large kinetic momenta /(q, 00) > /i(q, 00). Addi- 
tionally, the momentum spectrum is strongly enhanced around 
the maximum of /i(q, 00) (denoted as m\) due to the influence 
of the short pulse {£2, T2}. Moreover, the maximum of /(q, 00) 
(denoted as m) is shifted a little towards higher momenta com- 
pared to the pure Schwinger mechanism /i(q, 00), i.e. m > m\. 

In order to better understand this behavior, let us consider 
the influence of different fo values on the pair creation process 
in Fig. |4] For fo < (two upper figures), i.e. the maximum 
of the weak pulse occurs prior to the maximum of the strong 
pulse, we observe that the two curves coincide around m\. Due 
to the fact that the weak pulse acts prior to the strong pulse, 
this field configuration behaves much like a single strong pulse. 
The only difference is that particles which are created by the 
weak pulse are accelerated by the subsequent strong pulse such 
that m > mi. Oppositely for fo > (two lower figures), i.e. 
the maximum of the weak pulse occurs after the maximum of 
the strong pulse. In this case, the particles which are created 
by the strong pulse get an extra kick by the subsequent weak 
pulse such that f(q, 00) is shifted to larger momenta compared 



r =-T 2 /2 




Figure 4: Comparison of the distribution function /(q, 00) (solid line) with 
/l(q, 00) (dotted line). Depending on to, the maximum of the weak pulse occurs 
either prior or after the maximum of the strong pulse. The other parameters are 
chosen as in Fig. [3] 
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to /i(q, °o). Additionally, the pairs which are created by the 
weak pulse are accelerated less than in the previous case such 
that m < rri\. 

For to = we observe a combination of those two effects (cf. 
Fig. |3J: On the one hand we find that /(q, oo) > /i(q, co) for 
large momenta, i.e. the particles get an extra kick by the weak 
pulse. On the other hand, we also find that m > m\, i.e. the 
particles which are created by the weak pulse are still strongly 
accelerated by the strong pulse. 

Due to the fact that we obtain the maximal enhancement at 
to = 0, we consider this situation in the following. Before dis- 
cussing the enhancement of the number density n[e + e~] in the 
next section, we want to investigate the momentum spectrum 
depending on the pulse length T2 of the weak pulse in Fig. [5] 

The weak pulse becomes shorter by decreasing the value 
of T2, i.e., the energy of the corresponding quanta gets larger. 
Accordingly, it is getting easier to create an electron-positron 
pair via dynamical pair creation. As soon as the value of 
T2 < 5 • 10~ 6 <?V~', we observe distinct momentum signatures 
around the maximum m in tune with the inverse time scale 1 /T2. 
This might be interpreted as the onset of multi-photon domina- 
tion, i.e., whilst the dynamically assisted Schwinger mechanism 
was still non-perturbatively dominated in Fig. [3] the effect be- 
comes mostly perturbative in Fig. [5] 



m 

0.00020 - 



a[e*e-} 





Figure 5: Comparison of the distribution function /(q, oo) (solid line) with 
f\ (q, oo) (dotted line) for different values of tz. Top: T2 = 2 - 10~ 6 eV' 1 Bottom: 
T2 = 10~ 6 eV~ l The other parameters are chosen as in Fig. [5] 




Figure 6: Comparison of the number density n[e + e 
effects. The parameters are given in the text. 



according to the various 



3.2. Number Density 

We now turn to the investigation of the number density 
n[e + e~] for the dynamically assisted Schwinger mechanism. In 
Fig. [6] we vary the pulse length T2 of the weak pulse and, ac- 
cordingly, the combined Keldysh parameter y. The other pa- 
rameters are chosen to be {fii = 0.1£ c , Ti = 2 • 10 _4 eV -1 } 
for the strong pulse (this corresponds to j\ ~ 0.097) and 
{&2 = 0.01fi c , T2 variable) for the weak pulse: 

• The number density «j[e + e~] of a single strong pulse 
{£i , Ti } (pure Schwinger mechanism) is represented by the 
lower lying dashed-dotted line. 

• The number density «j +2 [e + e~] of a single strong pulse 
{fii + £2>T"i} (pure Schwinger mechanism) is represented 
by the higher lying dashed-dotted line. 

• The number density n™[e + e~~\ of a single weak pulse 
{£>2,T2} (pure multi-photon effect) with T2 being varied is 
represented by the dotted line. 



The number density n"" 2 [e + e ] which is obtained by 
naively adding the contributions of both the Schwinger and 



the multi -photon effect n s { [e + e ] 
by the dashed line. 



n [e + e ] is represented 



• The number density n^ +2 [e + e~] which is obtained by 
our numerical calculation of the dynamically assisted 
Schwinger mechanism is represented by the solid line. 

Additionally we have a dotted vertical line which corresponds 
to the parameter set at which the weak pulse has a time scale 
being equal to the Compton time of the electron, i.e., T2 = 1 /m. 
Note that these calculations have been done for q ± = 0, i.e. ig- 
noring the perpendicular momentum. Including the orthogonal 
momentum in the calculations is straightforward but not done 
for computational reasons. We observe the following behavior: 

• For a very small Keldysh parameter y <<c 1 such that 
T2 = Ti , we obtain the same result as the pure Schwinger 
mechanism with the parameters {£; + E>2,t{\. This is ex- 
actly what we expect since: 

E(t) = E x (t) + E 2 (t) = (fii + fi 2 ) sech 2 (f/n) (10) 
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Figure 7: Relative Enhancement of the number density n ll l+2 [e + e~]/n"" 2 [e + e~] 
due to the dynamically assisted Schwinger mechanism. The parameters are 
given in the text. 

• For a very large Keldysh parameter y » 1 the multiple 
time scale effect asymptotically coincides with the multi- 
photon result. This is reasonable because for very large y 
the multi-photon effect dominates the Schwinger mecha- 
nism by orders of magnitude. 

• In between these two extreme cases we have the region 
which is of special interest to us. Here, we see clear signa- 
tures of the dynamically assisted Schwinger mechanism: 
F° r 1 J ^ 2 we observe that n^" 2 [e + e~], i.e. the naively 
added Schwinger and multi-photon effect, is well below 
nf +2 [e + e~], i.e. the actual outcome of the dynamically as- 
sisted Schwinger mechanism. 

Actually, the relative enhancement of the number density 
nf 2 [e + e~]/n*"? 2 [e + e~] strongly depends on the various param- 
eters, as shown in Fig. [7] To this end, we vary &\ but keep the 
ratio between the field strengths £2 = 0.1£i constant. Addi- 
tionally, Ti = 2 • 10~ 4 eV _1 whereas T2 is varied again. 

The most important result is that in a region around y ~ 2 the 
relative enhancement reaches a maximum with an enhancement 
factor of up to 500 for a given parameter set. 

4. Conclusions and Outlook 

We have presented a detailed calculation for an example 
of the dynamically assisted Schwinger mechanism which has 
been suggested recently Si. This mechanism, i.e., electron- 
positron pair production from vacuum by a combination of laser 
pulses with different time scales, suggests that the rate of pro- 
duced pairs can be significantly enhanced by dynamical effects. 
First of all, we verified previous results, most prominently that 
one encounters a region around the Keldysh parameter y ~ 2 
where the relative enhancement reaches a maximum. Hereby, 
enhancement factors up to 500 have been obtained. Second, 
whereas the former investigations [fl |§| used the world-line 
instanton approach, the presented calculations are carried out 
within the framework of quantum kinetic theory. This method 
intrinsically enables us to additionally present valuable momen- 
tum space information of the generated positron spectrum. 



On the basis of these results we conclude that the search for 
the dynamically assisted Schwinger mechanism could be facil- 
itated by an ideal parameter choice in future experiments. It is 
obvious that the trivial combination of the two electric pulses 
considered here is neither optimal nor accomplishable experi- 
mentally. Therefore, a natural next step will be the investigation 
of more realistic electric pulses tuning the parameters for an in- 
creased pair production. Such pulse shaping in order to obtain 
an enhancement factor close to the greatest possible one even 
for more realistic electric field configurations is, although com- 
plicated, certainly feasible. Given the recent progress in theory 
and experiment a measurement of the Schwinger mechanism in 
this decade seems realistic. 
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